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Abstract 

The new approach to decays of heavy quarkonia, based on factor- 
ization and expansions in powers of the relative velocity of the valence 
quarks is reviewed. 



1 Introduction 

Heavy quarkonia are mesons, where both the valence quarks and the 
valence antiquark are heavy i.e. have masses much larger than Aqcd- 
Since the lifetime of the i-quark is too short for hadronization, only the 
discovery of the bc(cb) quarkonia is expected besides the well-known 
bb and cc quarkonia. Many features of the heavy quarkonia can be 
described by applying to the two-body system Q + Q the ordinary, 
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nonrelativistic Schrodinger equation with the interaction potential de- 
pending only on the relative coordinate r = \tq — Tq\. The earliest 
choice for this potential Q, which is still popular (cf. e.g. 0), is the 
Cornell potential 

V(r) = -- + br + c, (1) 
r 

where a, b, c are positive constants. Choosing the quark mass and 
solving the Schrodinger equation on finds the energy eigenvalues and 
the wave functions corresponding to the bound states. 

The energy eigenvalues are interpreted as quarkonium masses. The 
corresponding states are labelled like for the hydrogen atom IS, IP, 2S 
etc. Spin and the total angular momentum can be included in the 
usual way by writing 1 1 S'o, 1 3 Si etc. We shall not consider states, 
which are so heavy that they can decay strongly into pairs of heavy- 
light mesons containing the Q and the Q separately. Such quarkonia 
are broad and the simple description presented here becomes unreli- 
able. 

Using textbook formulae and the wave functions one can calcu- 
late the probabilities of the electromagnetic dipole electric and dipole 
magnetic transitions between various states of the quarkonia. Higher 
multipole transitions could also be calculated, but they have not yet 
been observed experimentally. Corrections to the usual formulae are 
sometimes introduced, because the ratios of the quarkonium sizes to 
the wave lengths of the emitted radiation are much larger than the 
corresponding ratios in atomic physics. Similar, but less reliable, cal- 
culations can be made for the n'S — * nS transitions (n < n') with the 
emission of a pair of pions. 

This approach has two difficulties, which one tries to overcome 
and a basic difficulty, which makes people look for completely new 
approaches. The first difficulty is that relativistic effects, e.g. fine and 
hyperfine splittings of the energy levels, are left out. This is overcome 
by using various "relativistic" equations instead of the Schrodinger 
equation and by introducing more general interaction potentials. Since 
no generally accepted two-body relativistic quantum mechanics exists, 
however, these approaches are strongly model dependent. They give 
the necessary splittings, but the centres of the multiplets are not re- 
produced any better than in the nonrelativistic theory. A comparison 
of 27 relativistic and nonrelativistic potential models with each other 
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and with experiment is given in the review article |J. A nonrelativistic 
potential, which reproduces the data for the bb quarkonia within their 
experimental errors and also gives good results for the cc quarkonia, 
has been described in Q]. 

Another difficulty is that quantum mechanics, strictly speaking, 
applies to processes, where the number and quality of particles do 
not change. A way out of this difficulty is to assume factorization. 
E.g. the probability of annihilation of an S-wave quarkonium into a 
lepton pair is evaluated as a product of two factors: the probability 
that Q and Q meet, this is just ^(O)) 2 , and the probability that then 
they annihilate into a lepton pair, which is calculated using quantum 
field theory. Such calculations are less reliable than the calculations 
of masses or non-annihilation transitions, but also they seem to give 
reasonable agreement with experiment. Incidentally, the analogy with 
calculations of decay constants of the heavy-light mesons (cf. |J and 
references given there) suggests that there may be difficulties here, 
which are hidden for the moment by other uncertainties. 

The main problem with potential models is, however, that their 
relation to sound theory is unclear. They also often produce incon- 
sistencies. E.g. calculating from the Schrodinger equation the ki- 
netic energy, one finds typically (v 2 ) rs 0.08c 2 for the bb(l 3 Si) and 
(v 2 ) 0.25c 2 for the cc(l 3 <Si). This suggest important relativis- 
tic corrections, which, however, are not contained in the Schrodinger 
equation. Attempts to calculate radiative corrections sometimes lead 
to infinities. One finds anomalous dimensions tending to infinity when 
~ — > i.e. in the nonrelativistic limit (for a discussion cf. |J) and 
non-cancelling infrared singularities in decay amplitudes J?]], ||, ||]. 

In the following we present a new approach Q , [|l0| , which explains 
at least some of the difficulties of the potential models and also brings 
a number of new interesting insights. 

2 Scales 

Let us consider the limit, when the masses of the valence quarks (i.e. 
of the valence quark and of the valence antiquark) are large and their 
velocities small. Here and in the following velocity means velocity 
in the rest frame of the heavy quarkonium. Thus, v <C c, or in the 
usual system of units, where the velocity of light c = 1, t) < 1. The 
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inequalities 



M » Mu » Mv 2 > . . . (2) 

where M is a mass of the order of the mass of the heavy quarkonium, 
or of the heavy quark, define a set of well separated energy scales. Let 
us discuss, which characteristics of the quarkonia are related to which 
scale @. 

The largest scale M corresponds to the total mass of the quarko- 
nium. It is also the inverse of the annihilation radius - the valence 
quark and the valence antiquark can effectively annihilate each other 
only when they are separated by a distance of order i or smaller. 
The proposal [lpj is to introduce into the theory a cut off A = M 
chosen so that for momenta below A the motion of the heavy quarks 
is non relativistic. The advantage of this choice is that heavy quark 
production, hard gluon emission and large relativistic effects in the 
motion of the heavy quarks are excluded. The price to pay is that 
new operators appear in the Lagrangian. This is like in the theory of 
electroweak processes, where for low energy processes a cut off exclud- 
ing the intermediate bosons can be introduced, but then four-fermion 
interactions must be added to the Lagrangian. 

The scale Mv is the scale of the momenta of the heavy quarks in 
the quarkonium. By the uncertainty principle its inverse is the scale 
of the quarkonium radius. The wave vectors of the gluons forming the 
average colour field within the quarkonium are also k ~ Mv. 

The scale Mv 2 is the scale of the kinetic energy of the heavy quarks. 
It is natural to assume that it is also the scale of the potential en- 
ergy and of the excitation energies of the quarkonia. It happens that 
h-QCD ~ 0.4 GeV is of the same order of magnitude. It is assumed 
that the stationary state vectors of the quarkonia consists not only 
of components \QQ), where the valence quarks are in a colour singlet 
state, but also of components \QQ,g), \QQ,gg), ■ ■ ■ , where the "dy- 
namic" gluons g have wave vectors k ~ Mv 2 . Note that in \QQ,g) 
the QQ system must be in the colour octet state in order to make the 
system QQg a colour singlet. 

Lower scales correspond to fine and hyperfine splittings, but we 
shall not discuss them here. 

Let us see now how these scales can be used to make various useful 
estimates. For very heavy valence quarks the radius of the quarko- 
nium is small and it is believed that the interaction is approximately 
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Coulombic. Since the kinetic energy and the potential energy should 
be of the same order, we find 



~ Mv 2 , (3) 

where R = jt- is of the order of the quarkonium radius. Eliminating R 
and using the fact that a s (fi) is a decreasing function of its argument 
we find 0] 



c ~ a s (Mv) > a s (M). (4) 

This implies that, at least in the high mass limit, it is inconsistent to 
include higher order corrections in a s (radiative corrections) without 
simultaneously including the higher order corrections in v (relativistic 
corrections) to at least the same order. 

In the Coulomb gauge one finds from standard perturbation theory 



\gA\ ~ a s (k)vk. (5) 

This shows that soft gluons couple weakly to slow quarks, even when 
a s is not small. It is remarkable that \gA\ ~ Mv s both for the gluons 
with k ~ Mv, for which a s (k) ~ v, and for the much softer dynamic 
gluons, for which k ~ Mv 2 , but a s (k) ~ 1. 

Let us estimate now the probability that within the quarkonium 
the valence quarks are in a colour octet state, which requires an addi- 
tional gluon to compensate the colour charge. It is assumed that this 
additional gluon is dynamic. We compare two rough estimates of the 
quarkonium energy shift due to the \QQ,g) component 



A g E~ P(QQ,g)AE, (6) 
A g E ~ (H\ [ d 3 xgv-A\H). (7) 



The first estimate, where AE ~ Mv 2 is an estimate of the energy 
of the additional gluon, is probabilistic. In the second, perturbative 
estimate, where H denotes the quarkonium, v-A ~ Mv 4 and the other 
factors do not change this estimate of A g E. Eliminating A g E we find 
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P(QQg) ~ v 2 . 



(8) 



Thus, contrary to the picture suggested by potential models, this prob- 
ability is small, but nonzero. 

Let us note one more estimate. Since the squared modulus of the 
wave function integrated over the volume occupied by the quarkonium 
(~ (Mv)~ 3 ) gives unity, ip ~ ¥2^2 in the region of interest. 



3 Lagrangian 

The full Lagrangian for the quarkonium is, of course, the ordinary 
Lagrangian of the standard model. Here, however, we need an effective 
Lagrangian with a cut off at A = M. The price for the cut off is 
that the Lagrangian has an infinite number of terms. In the small 
velocity limit, however, only a few of these terms survive. Systematic 
expansions in powers of velocity for the necessary matrix elements can 
be constructed by adding at each step only a finite number of terms 
to the previous approximation to the effective Lagrangian. In general 
the effective Lagrangian is written in the form 



C = Cught + Cheavy + SC. (9) 

Here 

Cu 9 u = -\ttG%G^ + Y,*jUhi (10) 

3 

is the Lagrangian for the gluons and for the light quarks assumed for 
simplicity massless. The second term 

Ueavy = (iD t + ^ * + X t (i Dt _ x (11 ) 

is the nonrelativistic approximation to the Lagrangian of the heavy 
quarks (created by and heavy antiquarks (created by x) interact- 
ing with the colour field. The mass M, which till now has been only 
esimated, is unambiguously defined by this equation. Everything else 
necessary to make the effective theory equivalent to the exact theory 
at the low energy scale is contained in the correction term SC. The 
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idea is to classify the infinitely many terms contributing to 5C accord- 
ing to their order with respect to velocity. Using information about 
the scales of various operators, as described in the previous section, 
one sees that Cheavy contains terms of order v 5 . Putting 5C = 0, one 
obtains an approximation, where in particular the numbers of heavy 
quarks and heavy antiquarks are separately conserved. This approx- 
imation may be used to estimate the masses of the quarkonia, but 
not their annihilations into light particles. The next approximation 
includes the bilinear operators of order v 7 



5 2 C Cl 



^ p L ^ t (^ 2 )V-xWx]+-", (12) 

where the dots replace the other terms with the same structure in the 
heavy fields. Three of them are of the same order v 7 . The coefficients 
Cj are calculable and made dimensionless by explicitly factoring out 
suitable powers of M. This approximation gives an improved descrip- 
tion of the structure of the quarkonia, but again does not allow the 
annihilation of heavy quarks. 

In order to describe annihilation one needs four-heavy-fermion op- 
erators, which occur in the next order in the expansion of 5C. One 
has 

4 

S 4 £ = M~ 2 f,> r ( ) j\\(>.i<- + ■■■, (13) 

j'=i 

where 0\ = 1, O2 = <?, O3 = A a , O4 = X a a. The fj are calculable, 
dimensionless coefficients and the dots denote higher dimension oper- 
ators with the same structure in heavy quark fields. The operators 
occurring in this formula are of order v 6 , but the coefficients fj are of 
order a 2 , therefore, S4C is of order v s . These terms in the Lagrangian 
make annihilation possible. Actually, the first term annihilates the 
1 S'o colour singlet states, the second the 3 5i colour singlet states, the 
third the colour octet 1 Sq states and the fourth the colour octet 3 S\ 
states. For annihilation of quarkonium H into light hadrons one finds 

F(H -» l.h.) = 2Im(H\5 4 C\H), (14) 

or approximately 
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4 

T(H - ZA) = 2M- 2 ^ /m/i (i^O^O^tf )• (15) 

3=1 

Note that the operators j = 3, 4 can annihilate only the quarkonium 
component 5), which has an amplitude of the order \J P(QQ, g) ~ 
v, and that that reduces the contribution to T by a factor of order v 2 . 
This usually justifies the omission of j = 3, 4 in low order calculations, 
but in the following section we show an example, where this reduction 
is compensated by another factor and the inclusion of j = 3,4 is 
crucial. 

For annihilations into leptons and/or photons the final state is a 
vacuum of hadrons, therefore, the operators j = 3, 4 cannot contribute 
and, moreover, e.g. 

{h\^ X xH\h) = \(h\^ x \q)\ 2 . (16) 

The same formula used for decays into light hadrons is only an ap- 
proximation known as the vacuum saturation. Rigorously 

{H\^xxH\H) = ]T \{H\^ X \X)\ 2 , (17) 
x 

where the vectors \X) are a complete set of states. Then the vacuum 
saturation means neglecting all the contributions with \X) 7^ |0). Let 
us estimate the terms, which are thus neglected. Terms with \X) = 
\g) do not contribute, because the operator ip^x does not change the 
colour and the state \H) is by assumption colourless. Terms with 
1-^0 = \99) contain only the components of \H), which besides QQ 
contain two dynamic gluons. Such components of H are of order 
v 2 , therefore, the corresponding terms on the right hand side of the 
previous formula are reduced by a factor of order 0(v A ). For \X) = 
\qq) or \X) = \ggg) the reduction is even stronger. Thus, in the present 
approach the error inherent in the vacuum saturation approximation 
can be rigorously estimated and is small. 

4 Examples 

As our first example let us consider the annihilation of rj c into light 
hadrons. The leading term is 
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r(r? c -> ZA) 



M 2 



(18) 



Using the vacuum saturation approximation as explained above 



The matrix element on the right-hand side can be related to a radial 
wave function averaged over a volume of order M~ 3 



where the factor three under the square root is the number of colours. 
Using again the scale estimates one finds 



Thus, up to corrections of order v 2 the decays of corresponding 1 Sq 
and 3 Si states are governed by the same matrix element. This is the 
realization in the present approach of the spin symmetry from the 
heavy quark effective theory. At the same level of precision one can 
connect the matrix elements for the production and for the decays of 
the heavy quarkonia fl(j[| . 

Let us note an advantage of the present approach as compared with 
potential models. Certain expressions, which are infinite in potential 
models, like V 2 R(r) for r — > 0, can be here regularized by standard 
methods of quantum field theory. 

As our second example let us consider the decay of a P-wave 
quarkonium into light hadrons. In potential models this problem was 
plagued by infinities. We shall show how it is resolved here. Let us 
note first that for P-states ip(0) = 0. Therefore, choosing for defi- 
niteness the 1 l P\ quarkonium h c as our example and identifying the 
wave function Rh c with the wave function given by potential models 
we have 




(19) 




(20) 



R 1 p(0) = R Vc (m+O(v 2 )). 






The next term contains 



{h c \i>\-\D)xx\{-\D)i>\h c ) 




(23) 
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where R' hc is the derivative with respect to r of the radial wave function 
of the valence quarks in the quarkonium h c . The right-hand side is 
well known from potential models. This matrix element, however, is 
suppressed by two powers of v as compared to the matrix element 
responsible for the decays of S-w&ve quarkonia into light hadrons. At 
the same order in v there is another matrix element 

(h c \^X a xx^ a ^\h c ), (24) 

which corresponds to the annihilation of the \QQ,g) component of h c 
and which has no analogue in potential models. The suppression by 
v 2 due to the choice of the \QQ,g) component of the h c state vector 
is compensated by removing the factor of order v 2 due in the matrix 
element ( |23| ) to the two factors D- To put it more intuitively, the 
gain due to the fact that the colour octet QQ system is in an S-state 
and the quark and antiquark can easily meet to annihilate compen- 
sates the penalty for introducing the additional dynamic gluon. When 
both terms ( p3[ ) and ( |24| ) are included, one finds that the infrared di- 
vergencies of the radiative corrections cancel and thus a meaningful 
calculation can be performed. One can also use the renormalization 
group equations to improve the results fllPf . 

5 Conclusions 

The problem of heavy quarkonia is difficult. Brute force expansions 
in a s or in i are unlikely to work. Simple-minded potential models 
work well within their applicability range, but the reason for their 
success is not understood and obvious contradictions are encountered 
on the way. Also the parameters of these models (e.g. quark masses) 
have no definitions sufficiently precise to connect them to quantities 
occurring in QCD (in the case of quark masses to MS or pole masses). 

A recently proposed combination of factorization and expansions 
of parts of the matrix elements in powers of the relative velocity v 
H [0] seems very promising. It is much closer to rigorous QCD than 
the potential models. Its various assumptions and approximations can 
be rigorously studied, because the relevant quantities are defined in a 
way, which is understandable from the point of view of QCD. Its weak- 
ness at present is that it introduces many matrix elements, which in 
principle are calculable from lattice QCD, but in practice may remain 
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for many years to come just phenomenological parameters. It should 
be stressed, however, that the theory already has given some interest- 
ing approximate relations between these matrix elements and that it 
removes the inconsistencies, which have been plaguing the potential 
models. 
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